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We consider the asymmetric simple exclusion processes (ASEP) on a ring con- 
strained to produce an atypically large flux, or an extreme activity. Using quantum 
free fermion techniques we flnd the time-dependent conditional transition probabil- 
ities and the exact dynamical structure factor under such conditioned dynamics. In 
the thermodynamic limit we obtain the explicit scaling form. This gives a direct 
proof that the dynamical exponent in the extreme current regime is z = 1 rather 
than the KPZ exponent z = 3/2 which characterizes the ASEP in the regime of typ- 
ical currents. Some of our results extend to the activity in the partially asymmetric 
simple exclusion process, including the symmetric case. 

I. INTRODUCTION 

The Asymmetric Simple Exclusion Process (ASEP) is a minimal model for many traffic 
and queueing processes [l|-l3| . It describes a non-equilibrium system of many driven particles, 
interacting via hard core repulsion. This Markov process is defined on a lattice, each site 
of which can be empty or occupied by one particle. Particles jump independently after an 
exponentially distributed random time with mean ^/{p + q) to a nearest neighbor site on 
the right or on the left, provided that the target site is empty (hard core exclusion rule). 
The probability to choose the right neighbour is p/{p + q) while the probability of choosing 
the left neighbour is g/(p + g). In spite of its simple formulation, and very simple product 
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stationary state, the time- dependent characteristics are very nontrivial to obtain, even in the 
simplest periodic one-dimensional case (closed ring of L sites). This is due to the fact that 
e.g. the time- dependent operators of type e^*/(0)e~^* are complicated objects because the 
generator of the stochastic dynamics H is equivalent to an interacting many-body quantum 
Hamiltonian. 

Recently, it has been found that the ASEP, under the restriction to produce an extreme 
flux or extreme activity is governed by effective long range interactions ^|. In this setting 
one considers realizations of the process for a duration T which for a long interval between 
large times t and T — t (where T — t is itself also large) have carried an atypically large flux. 
This extreme event quite surprisingly makes the conditioned process intrinsically related to a 
much simpler system than the original one, namely to a system of non-interacting fermions. 
This fact suggests to use field theoretic free fermion techniques to compute time- dependent 
correlations in this particular large deviation limit of the classical stochastic dynamics of 
the conditioned ASEP. 

In the present contribution we shall employ these techniques to calculate two particular 
time-dependent correlation functions: the conditional probabilities P(?7, t; ?7o, 0) for the tran- 
sition from a microscopic A^-particle initial configuration rjQ to a final configuration r] at time 
t and the dynamical structure factor which is the Fourier transform of the time- dependent 
density-density correlation function (?T-fc(t)?T-o(0)) — in the stationary ASEP in the two 
above mentioned limits: under the restriction to produce a very large flux and under the 
restriction to produce a very large activity. We note that for the totally asymmetric version 
of the ASEP (TASEP) on an infinite lattice and on a ring the conditional probabilities were 

n n 

found in [5| and [6| respectively using the Bethe ansatz. The dynamical structure for the 
infinite system was obtained by Prahofer and Spohn using random matrix techniques 3]- 
We remark that that the time- dependent correlations for atypical stochastic trajectories (in 
the large deviation limit) are practically unaccessible by numerical methods, even though 
some promising approaches have been developed recently Q]. 

The paper is organized as follows. In Sec. [Ill we introduce the ASEP conditioned to carry 
a certain nontypical flux or to exhibit a nontypical activity. In the first two subsections we 
review for self-containedness the our previous work ^ which among other things clarifies 
how conditioning on large flux and activity is related to free fermions. In the third subsection 
of Sec. [TTlwe present a new result, viz. the form of temporal correlations in the large current 
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or large activity regime. In the following main sections we focus on large flux and derive the 
conditional transition probabilities (Sec. IIII|) and dynamical structure factor (Sec. II Vp . We 
conclude with a brief summary and some remarks on the large limit activity in the ASEP 
and symmetric exclusion. 



II. ASEP CONDITIONED ON A LARGE FLUX 



A. Definitions and relation to an exclusion process with long-range interaction 

The flux (or synonymously the current) in the ASEP is the mean net number of jumps 
in some time interval. Conditioning the ASEP to carry an typical flux is achieved by first 
ascribing to the ASEP space-time trajectories weighting factors e'^"^ where s is a generalized 
chemical potential, conjugated to the total current J, registered along the trajectory up to 
time t (see e.g. jo] for a elaborate treatment). Fixing s, one picks up an ensemble of space- 
time trajectories characterized by any desired average current j{s). In particular, fixing the 
desired current to be very large (in the sense described in the introduction), one obtains a 
particular set of trajectories, where several tendencies are clearly seen {4] when taking the 
limit s — )■ 00: (a) the backward hoppings are totally suppressed (b) formation of clusters 
of particles, potentially leading to the flux reduction, is strongly suppressed (c) long range 
interactions between the particles are generated. 

In more details, it was noted that the ASEP under the restriction to produce extreme 
flux after proper rescaling of time is described by a Master equation 

d\P) 
dt 

with effective stochastic Hamiltonian 



-Heff\P) (1) 



iJe// = Ai/A-i-/i, (2) 

where A is a diagonal matrix with positive entries given further below, and n is the lowest 
eigenvalue of the Hamiltonian H 

L 

H = -Y.a^a;,^,. (3) 

n=l 

with the spin 1/2 lowering and raising operators of the Lie algebra SU (2). We point out that 
in the limit of very large flux the same stochastic process describes the partially asymmetric 
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simple exclusion process (ASEP) or even the symmetric case, since in this limit any backward 
hopping does not contribute. 

Similar considerations can be carried out for the activity where space-time trajectories 
are weighted by a factor e'^'^ where A is the total number of jumps along the trajectory, 
irrespective of the direction of the jump. In the case of ASEP restricted to have an extremely 
large activity (counting hoppings, irrespective of the direction), the effective Hamiltonian 
has the form ([2]) where H is substituted by 

L 

Hn = -Yl (P^n f^n+l + Q(^n+l(^n) . (4) 
n=l 

and /i is substituted by the lowest eigenvalue of Hjj and p is the hopping rate to the right 
(clockwise) while q is the hopping rate to the left (anti-clockwise) of the ASEP. Both Hamil- 
tonians ([3]),(j4]) are free fermion Hamiltonians, and in particular for the case of symmetric 
hopping p = q = 1, the so called Symmetric Simple Exclusion Process, or SSEP, the dynam- 
ics is governed by the effective Hamiltonian ([2]) with H substituted by the physical XXO 
Hamiltonian describing quantum fermions on a ring of L sites, 

L 

HxXO = -J2 i'^n^n+l + ^n+l^n ) (5) 
n=l 

and fi substituted by the lowest eigenvalue of Hxxo- All the above mentioned Hamiltonians 
commute with each other and share the same ground state The components of in the 
natural basis constitute the respective entries of the diagonal matrix A from ([2]) (rjlfi) = A^^. 
By the Perron-Frobenius theorem the ground state is nondegenerate, all its components are 
positive (?7|yu), so the matrix A~^ exists. It is straightforward to verify that the stationary 
state (the state with zero eigenvalue) of the stochastic process ([I]) is given by A|/i). The 
stationary state is thus common for all Hamiltonians mentioned above, for details see |4l|. 10 1. 

It may be worthwhile noting that the effective process is an interesting object of study 
in its own right, i.e., without reference to atypically large currents in the usual ASEP. This 
process is a totally asymmetric exclusion process with long-range interaction. It describes 
nearest neighbour hopping, but the hopping rate depends on the positions of the other 
particles. Specifically, the move of the k-th particle, located at position Uk in a configuration 
rj, to the consecutive position + 1 in the new configuration rj' has the rate jj] 

f^^ sm(7r(rafc - ni)/L} 
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where the product is over / and k from 1 through A^. If the target site + 1 is occupied, 
the hopping rate is zero, in agreement with the exclusion rule. 

For the case of extreme activity one obtains in the same fashion a partially asymmetric 
or symmetric exclusion process with long range interaction. The symmetric case has been 
studied some time ago by Spohn 15[ who derived a hydrodynamic fluctuation theory for 
the large scale dynamics of the process. Noticing that fluctuations are driven by a Gaussian 
field he presented a general form of the dynamical structure function which represents a 
universality different from the Edwards Wilkinson equation. He also noted the link of the 
hopping rates (Ej) to Dyson's Brownian motion, which was pointed out independently in ^ 
for the totally asymmetric case. Thus our discrete model has a natural interpretation as 
a discrete, biased random walk version of Dyson's Brownian motion driven by an external 
field. 



B. Free Fermion states and spectrum 

The complete set of eigenstates of H, Hjj and Hxxo in a sector with particles is 
characterized by a combination of plane waves where each plane has a quasimomentum a 



Therefore each 



satisfying e^°'^ = (—1)^"''^ because of the periodic boundary conditions 
a takes one of the quantized values 

ak = —ik —), A; G {!,..., L} (7) 

and thus each eigenstate is defined by an A^-tuple of integers {k} = {ki, . . . , kjy} where each 
kj G {1, . . . , L}. With these definitions we can write the eigenvectors as 

= ^ Xmim2...mNi{k})\mim2....mN) (8) 
mim2...»njv 

where the vector |mim2....m7v) denotes a state with particles at respective positions rrij (not 
necessarily ordered) and the sum is over all possible particle positions rrij G {1,2,...,L}. 
The free fermion wave function in coordinate representation has the form 

Q 

where Q is a permutation of indexes 2, A^) = {Qi,Q2, ■■■iQn) and (—1)*^ denotes 
the sign of the permutation. The factor T{m} is zero if in the set {mim2....mAr} some rrij 
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coincide, and otherwise is equal to 1 or —1, depending on whether the set {77111712 ■■..mjsf} 
was obtained from the ordered set by an even or odd number of pair permutations. Any 
exchange of particle positions or quasimomenta results in the eigenfunction IfJ^ik}) changing 
sign, which reflects the Pauli principle. The states with all possible sets {k} are 

normalized, orthogonal and form an orthonormal basis in the respective sector of Hilbert 
space with number of particles. The eigenvalues E^''^ corresponding to the are 
sums of those for each quasiparticle 

N 

^w = E^^K) (10) 

where the quasiparticle energies are 

e\ak) = -e-"'" (11) 
£^^(«,,) = -(pe-'°'=+re'"'=) (12) 
£''''°(afc) = -2cosafc (13) 

for the Hamiltonians (I3]),(!l]) and ([5]), respectively. 

The ground state |yUo) which minimizes the energies E^f^y is characterized by the choice 
of quasimomenta where kj = j, i.e. for the set 

{k}o = {l,2,...N}. (14) 

For this choice the ground state energies Eq for all three cases are real. In the following 
by default we treat the Hamiltonian (j3]) with quasiparticle "energies" f lTT]) and drop the 
superscript /. 



C. Correlation Functions 



Consider a general time-dependent correlation function {G{t)F{0))eff in the stationary 
distribution of the effective stochastic process defined by ([2]). Here G and F are functions 
of the occupation numbers and hence represented in the quantum Hamiltonian formalism 
by diagonal operators and G{t) = exp{Hefft)G exp{—Hefft). The basis of our present work 
is the following simple, but fundamental property of conditioned processes of the form ([21). 

Theorem. Let the operators G,F be diagonal in the natural basis defined by i^e//- 
Then 

(G(t)F(0))e// = mt)F\^,) (15) 
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where G{t) = exp{Ht)G exp{—Ht) and H is the respective free fermion Hamiltonian ([3]),(jl]) 
or 

Proof. The proof is simple, but for readers not famihar with the machinery we provide 
here the details. By definition, 

(G(t)F(0))e// = (s|e^^//*Ge-^-^^*F|P*) (16) 

where the summation vector {s\ and stationary distribution vector \P*) are left and right 



lowest eigenvectors of -f^e// with eigenvalue 0. We recall that \P*) = A|/i) j4| and that 
the summation vector {s\ = (1111. is a vector with all unit components. The property 
{s\Heff = is simply a stochasticity condition (conservation of probability) for the Hamil- 
tonian Heff. Using the latter condition and the fact that [A,F] = when both A and F 
diagonal operators we obtain 

(G(t)F(0))e// = (.|Ge-^^//*AF|/i). (17) 

Now we note that e~^'=ff^A = Ae~^^~^^^, which is verified by expanding the exponent on 
both sides in Taylor series and comparing the series term by term, using ([2]). Inserting this 
equality into (fTTl) . and using [A, G] = 0, we obtain 

(G(t)F(0))e// = (.|AGe-(^--)*F|/.). (18) 

Finally, noting (s|A = and (/i|e^* = (/i|e^*, we obtain the right hand side of (fTSjl . □ 
Some remarks are in order. Theorem (fT5l) is straighforwardly generalized to 
multipoint space-time correlation functions and one obtains {G(ti)F(t2)-..Q(tn))eff = 
{lJ,\G{ti)F{t2) . . . Q(tn)\fJ'), provided that all operators G,F,...,Q are diagonal in the nat- 
ural basis, i.e., represent observables of the classical stochastic process generated by Heff. 
We also point out that the theorem is phrased for our specific case of the current or activity 
in the ASEP. However, a similar results applies for any conditioned stochastic dynamics in 
which case H = H{s) is the weighted generator of the unconditioned process. 

III. CONDITIONAL PROBABILITIES 

We pick two arbitrary configurations ?7o and rj of our system containing particles at the 
positions rii and mt respectively, i.e., |?7o) = \nin2.-..nN) , \ri) = \mim2....miq) . We assume 



8 



both sets to be ordered, < rij+i and rrii < rrii^i. The probabihty to find the system in 
configuration t] at time t, provided it has been in configuration 770 at time t = is given by 

ftfe<;*.o) = ^^^. (19) 

{V0)eff 

where the subscript denotes that the average is computed with respect to the stationary 
state defined by the effective dynamics ([1]). The operators 17(0) = |?7)(?7| and 170 = |?7o)('7o| 
are diagonal operators, which allows to use the Theorem (IT^ . The denominator in (1191) can 
be obtained by using f[T^ with G = I and F = fjo. Using the results for the stationary 
probabilities of ref. we find 

2N{N-i) n - n 

{Vo)eff = {M{Vo\fi) = mXoXvo = ^jr- n sin^TT^-^ (20) 

i,j 

l<i<j<N 

The right hand side of fl20|) are (modulo squares) the components of a Slater determinant 
with quasimomenta ( !T4|) chosen so as to fill the Fermi sea. For details, see e.g. 
For the numerator we have, using ([T5l) . 
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12|. 



(^(i)^o)e// = (/ile^V'^'^ol/i). (21) 

Since the eigenvectors ([8]) form an orthonormal basis the sum of projectors 
(A^!)"-*^ = / over all A^-tuples {k} is a unit operator. The factor (A^!)~^ 

appears because in the sum ^{fc} /({ctfc}) = X]j'i=i ••• X]^jv=i ^('^ii' '^h^ ■■■^jn) ^^ch different 
set of as occurs N\ times. Inserting it in (12T1) . we obtain 

^(/x|e^V-''*Elm)(ml^o|/x) = ^J2e^^'^-'^i''y^'Z*{r^)Z{vo), (22) 
W ■ {fc} 

where Eq is the ground state energy corresponding to the specific set of quasimomenta 

{k}o given by ([H]), and Z{ri) = {i2{k}\v\l^) = {N\yx*^{{k})Xr,{{k}o) are found using 1^. 

Substituting Z^t]) into (!22l) . one obtains 

m)Vo)eff = (iV!)'x;({A:}o)x.o({A^}o)e^*5^e-^{^>*x;({fc})x.({^}) (23) 

{k} 

Using the explicit form of x ©, the sum over {k} in fl23l) can be rewritten as 

^^'^ {k} Q,Q' 
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In the sum over the permutations ^qq/ there are (A^!)^ terms; however, under the sum- 
mation over ki, k2, ■■■k]^ and reshuffling of k-s only A^! terms are independent, namely. 

Substituting this into ( l24l) . we get 

L_ ^(_l)Qe^Ef=i(™c3,-n,)a,^-iEf=ieK^,)_ ^25) 
{k} Q 

So far the discussion has been general and applicable to all three Hamiltonians ([3]),(11D and 
([5]). Focussing now on the default case ([3]) we use f lTOj) and ( ITT]) for further simplification. We 
expand the exponents e"*^'-"*^ = ^^q^*^*"''*/'^' (125!) . and collect the terms with the same 
Ok- Under the summation over {k}, only the terms with niQ^ — rij — s = —kL contribute to 
the sum, where k = 0, 1, 2, ... if mq. > Uj and k = 1, 2, .. otherwise. Since each satisfies 
giOfc-L _ (^_Y^N+i ^ q^i^ number of particles N each contributing term in fl2S]) is equal to 
^1^1 1 = . For even iV, the signs of the terms with different k will alternate. Now we 
introduce the function 

9L{d, t) = J2 [{-irszgn{d)f^' (26) 

where d is an integer ranging from —L + 1 to L — 1 and di = d for d > and di = d + L for 
d < 0. The function sign{d) = 1 for d > and signed) = —1 for d < 0. With this function 
the expression ( 125|) can be rewritten in compact determinantal form as 

-j;^Y.^-^^'^U9Li^Q^-^k,t) = det[gL{mj-n,,t)] (27) 

Q k=l ^ 

where we used Leibnitz formula for the determinant. Finally, using (l20l) . ( 123|) and ( 1271) . the 
conditional probabilities f|T9l) can be brought after some algebra into the final form 



Pl{v. t; Vo, 0) = e^«* W det[gLimj - m, t)] (28) 

V WO/e// 

The (?7o)e//7 {v)eff are stationary probabilities of the initial and the final state respectively, 
given by f l2Q|) . At the last step of the calculation we have used the fact that all the compo- 
nents Xrji{^}o) of the ground state eigenvector can be made positive, see the discussion after 
Eq. (ED. Consequently x;(({/';}o)/x;o(({^}o) = Xvii{k}o)/Xmii{k}o) = V (v) ef f / {r]o) ef f ■ 
Eq. (l28l) is the main result of this section. The determinantal structure of the conditional 
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probabilities was first noticed by Spohn 15| for the symmetric case. In contrast to the sym- 
metric case, here the matrix elements of the determinant are the propagators of a totally 
asymmetric random walk. Extending the link of the symmetric model to Dyson's Brownian 
motion we may regard our model as a totally asymmetric Dyson random walk. 

Let us discuss some limiting cases of ( l28l) . For t = one has det[gL{dij,t)] = 5nr,Q, yielding 
the correct normalization P(?7, 0; ?7o, 0) = In the simplest case of one particle = 1, 
Eq = -1, (?7)e// = (?7o)e// = 1/-^, and we obtain 

PL{m, t- n, 0) = e-'gL{d, t) = ^ (29) 

where d = m — n for m > n and d = L — n + m otherwise, and n, m are initial and final 
particle positions. Each term e~^t'^~^'^^ / {d + kL)\ in the sum is a contribution of a Poisson 
process e~^\^ / k\ where an event (hopping of a particle to the right with rate 1) has happened 
d + hiL times during time t. Along this trajectory a particle starts from site n, and arrives 
at site m after making k complete circles on the ring of size L. Thus the parameter k in 
(126|1 is the winding number. 

The formula fl28|) has been obtained for the continuous time Markov process ([1]). For the 
discrete time update [13[ we expect the Poisson distribution terms in ( l26l) to be substituted 
with the binomial distribution terms 

^DISCRETE^^^ t) = Y,^ [{-lYsign{d)f^^ (^^ _^ ^ J _ py-dL-^L^ ^3Q) 

where p is the probability of a TASEP particle to hop. Unlike the function fl26|) . the sum 
( 130|) is truncated for any finite t. 



IV. DYNAMIC STRUCTURE FACTOR 



The dynamic structure factor in the large current regime of a periodic chain with L sites 
is defined as the Fourier transform of the stationary correlation function hiini, ti]n2,t2] p) — 

= {nn^{ti)hn^{t2)) ef f — P^, whcre fikif) = e^^ff^nke~^''-ff* are particle number operators 
and p = N/L is the stationary particle density. Without losing generality, we assume 
n2> rti. Because of the translational symmetry and time independence of the Hamiltonian 
Heff, the correlation function h depends only on the differences 77,2 —rii = n, ^2 — ^1 = ^, i-e. 
/i(ni, 77.2, ti, ^2) = hL{n, p, t). Thus we can write the real-space representation of the dynamic 
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structure factor as Siin.p.t) = hiin.p.t) — p^. Moreover, by particle-hole symmetry we 
have that Siiji, 1 — p,t) = SL{—n,p,t) and trivially 5*^(72, 0,t) = Si{n, l,t) = 0. Hence we 
can limit our discussion to the range < p < 1/2. In order to simplify notation we drop 
the dependence on p in the structure function and the correlation function and simply write 
SL{n,t) and hL{n,t). 

The number operator is diagonal and therefore the formula (fT5!) is applicable. hL{n,t) 
can then be calculated in similar manner as the conditional probabilities in the previous 
section. We shall present only the final result (the derivation proceeds analogously to the 
respective quantum mechanical calculation of (o"^2(^)^ni(0)) |l2|). 

iV L 1 ^ ^ 
^) = __ g-^afc'^+£(afe)t g^a^'^-£(a^ )* Qia^n-e(ak)t ^-ia,n+e{ai)t ^g-j^^ 

^ k=l 1=1 ^ k=l 1=1 

where the are of the form (I7j) with the ground state choice kj = j, and e{ak) is given by 
one of expressions ( fTTIfTSl) . 

Notice that even though the expression flM]) formally looks like the corresponding quan- 
tum formula in imaginary time, its analytical properties and limits are crucially different. 
All the sums in ( 13T]) are strictly real which can be verified straighforwardly, using sym- 
metricity of the ground state set f|T^ of pseudomomenta a^. We remark that for our default 
choice e{ak) = — e"*"*, the summation from 1 through L in Eq. lpTl) attains a simple form 
by expanding the exponent. We obtain then 

^5^e^"'"-^("')*=^7i(n,t), (32) 
1=1 

where gL{n,t) is given by fl26|) . 

From the two-point correlation function fl3T]) we compute the dynamic structure factor 

L 

SL{v.t) = Y,e-^^''^''"'SL{n,t) (33) 

n=l 

with the integer momentum variable p G {1,2,..., L}. Obviously Sl{Pi t) = Sl{p + nL, t) 
for any integer n and 5*^(0, t) = which allows us to restrict the subsequent study of the 
dynamic structure factor to the range pG{l,2,...,L — 1}. To evaluate fl5^ in this range 
we first observe that the Fourier transformation turns the exponentials of the summation 
variables a^, O-i into the Kronecker-delta 5p^k~i ■ Then we write the second sum in (l3T|l 
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(which runs up to A^) as a sum from 1 to L and subtract the part from + 1 to L. This 
yields as an intermediate expression 

1 ^ 



k=l 

N L 



p,k~l 



k=l l=N+l 



(34) 
(35) 



■.= s\:>{p,t)+s['>{p,t) 

for which we analyse next the double sum S)^ {p,t). We focus on the default case (ITTil . 

For the default case the difference of relaxation times in the exponential takes the simple 
form 

e{ak) - e{ai) = -(1 - e^^'P/'^)e{ak) (36) 
and for notational convenience we introduce 



tp := (1 - e^'''P/^)t. 



(37) 



Bearing in mind the range of definition of the momentum variable p G {1, 2, . . . , L — 1}, the 
Kronecker-delta in conjunction with the summation limits of the double sum gives rise to 
three distinct regimes for p. Careful analysis yields 

p 

L 



1 ^ 



k=l 

N 



1 



L 

1 
L 



p=l,...,N-l 
p = N,...,L- N 



(38) 



k=l 



N 



gtpe -fe p = L- N +1,...,L-1. 



y k=N+i-L+p 

In the thermodynamic limit L ^ 1 the sums over k and j in f l3ip turn into integrals, 
R{n,p,t) := limL^ooxEfcLie^""""^^"'^^* = ^ J^^^e'P''~<P^^dp. Notice that here p is real- 
valued and we define it to be in the interval [— tt, tt]. This yields 

S{n, t) := lim ^^(n, t) = R{n, 1, t)R{-n, p, -t) - R{n, p, t)R{-n, p, -t). (39) 

where i?(n, 1, t) = t"'/n\ corresponds to the limit of the function g{n, t), where only the first 
term in (!26|) appears (no winding condition). With a view in large, but still finite systems we 
point out that this observation imposes obvious validity limitations on the integral expression 
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(a) (b) 

Figure 1: Dynamic structure factors (ra.m(t)no(0)) — for p = 0.5 as function of time for < ?i < 8, 
computed from Eq. (j39p . Panel (a): connected correlation function h{n, t)—p'^ for odd n = 0, 2, 4, 6, 
represented by bold, thin, dashed, and dotted lines respectively. Panel (b): connected correlation 
function t) — p^ for odd n = 1, 3, 5, 7( bold, thin, dashed, and dotted lines respectively). 



( 13911 as an approximation for large system size: The limiting behaviour cannot be used as 
an approximation for t>L and L. The integrals R{n,p,t), apart from obvious special 
cases R{n,l,t) = t"'/n\, and R{n,p,0) = (nvr)"^ sin nvrp are not expressed in elementary 
functions and must be evaluated numerically. In Fig. [1] we show the function S{n, t) for 
even and odd n for the particular case of half- filling N/L = 1/2 at different times t. For 
t = the difference in the expression Sq for energies the quasiparticles become irrelevant and, 
using R{n,p,0) = (nvr)"^ sinnvrp, we obtain from ( l39l) the static density-density correlation 
function 

sin^nvrp , 

first derived in jl4l |. 

In order to explore the large-scale behaviour of the dynamic structure factor we study the 
behaviour for small momentum p and large times t. To this end we return to the definition 
f l55]l and first observe that in the thermodynamic limit 



where now tr, 



e*^)t. Using; t_ 



pn 



dx 



(41) 



-e *^tp allows us to rewrite this expression as a 



difference of integrals over the same function exp(e "tp) with integration intervals [— pvr + 
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p, pir + p] (positive term) and [— pvr, pn] (negative term) respectively. On the other hand 









dx 


2^ J 


f-fm 






' dx e*"*" 


2^ J 


— pTV 

r-jy-K 




' dx e*" 


. 2^ J 





p E [0, 2p7i] 

p e [-TT, . . . , -2p7r] U [2p7r, . . . , tt] (42) 
^-^^ p G [-2p7r, 0] 

Putting everything together we finally obtain the dynamic structure factor for p < 1/2 

fpTT+P 



1 

2^ 

Sip,t)={ i- 
Ztt 

1 

{ 27r 



— pTT 



rfxe*''^"" p G [0,2p7r] 

e*p''~'" p G [-vr, . . . , -2p7r] U [2p7r, . . . , tt] 
dxe^^"'" p G [-2p7r,0] 



(43) 



This, along with the symmetry relation S{1 — p,p, t) = S{p, —p, t) provides an exact integral 
presentation valid for all densities p, momenta p and times t. The static structure factor 
takes the simple form (p < 1/2) 



, — n G h2p7r, 2p7rl 
p G [-vr, . . 



(44) 



p pE [-vr, . . . , -2p7r] U [2p7r, 

cf. the real-space result pUj) . 

We are particularly interested in the large scale behaviour as expressed in the scaling limit 
of small p and large t of the form p^t = u where z is the dynamical exponent and u is the 
scaling variable. In the limit p — )■ only the first and third expression in f l43p are relevant. 
From the occurrence of the factor tp = {1 — e*^)t we conclude that there is non-trivial scaling 
behaviour for 2 = 1, i.e. for u = pt. In this scaling we have tp = —iut which yields the 
desired result 



S{u) 



\u\ 
27rt^ 



-iu cos piT—\u \ sin pir 



(45) 



which is valid for all p G [0 , 1] and in agreement with the universal form the dynamic 
structure factor derived in [l5| for the symmetric case. The presence of the particle drift does 
not change the universality class as it does for the usual unconditioned exclusion process 
where the undriven model is in the universality class of the Edwards- Wilkinson equation 
with dynamical exponent z = 2 while the driven model is in the KPZ universality class with 
z = 3/2. This is in agreement with an earlier observation that for stochastic dynamics which 
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have an underl ying; free-fermion structure an external drift can be absorbed into a Galilei 



transformation 
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V. FINAL REMARKS 



We obtained analytically conditional probabilities and the two point time-dependent den- 
sity correlation functions for the ASEP conditioned to carry a very large average current. 
The conditional probabilities have determinantal form and can be expressed through ele- 
mentary functions. The density correlation functions are obtained both for a finite system 
and in the thermodynamic limit. By Fourier transformation we have computed the exact 
dynamical structure factor and derived its large-scale behaviour. The natural scaling vari- 
able turns out to be u = pt which proves that the dynamical exponent of the conditioned 
ASEP is z = 1. From the explicit scaling form we read off the collective velocity Vc = cos pir 
of density fluctuations which is in contrast to Vc = {p — q){l — 2p) of the usual ASEP in the 
regime of typical currents. The relaxation part is symmetric in u and very different from 
the corresponding quantity in the usual ASEP Q which has dynamical exponent 2 = 3/2 
for the driven case and z = 2 for the symmetric case. In our model the presence of a drift 
does not change the universality class. 

Our results have a natural generalization to the study of large activity, i.e. to choosing the 
Hamiltonians @ and One has to replace the energies e{ak) in ( ET\i by the respective 
expressions — (pe"*"*" -l-re*"*^) and — 2cosa;fc. It will be interesting to study not only the 
hydrodynamic limit, but also the microscopic structure of shocks for the general case. It 
would also be interesting to extend the analysis of the effective dynamics of driven systems 
under large deviation constraints to the non-diagonal case, e.g. to compute current-current 
time-dependent averages. 
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